Transconductance fluctuations as a probe for interaction induced quantum Hall states 

in graphene 
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Transport measurements normally provide a macroscopic, averaged view of the sample, so that 
disorder prevents the observation of fragile interaction induced states. Here, we demonstrate that 
transconductance fluctuations in a graphene field effect transistor reflect charge localization phenom- 
ena on the nanometer scale due to the formation of a dot network which forms near incompressible 
quantum states. These fluctuations give access to fragile broken-symmetry and fractional quantum 
Hall states even though these states remain hidden in conventional magnetotransport quantities. 
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Progress in graphene sample quality starts to disclose 
rich physics related to Coulomb interactions as well as 
interaction induced lifting of symmetries associated with 
the spin and pseudospin degrees of freedom. This is vis- 
ible in macroscopic magnetotransport studies as addi- 
tional incompressible quantum Hall condensates. If inter- 
actions are ignored, quantum Hall plateaus are expected 
and seen in experiment for the filling factor sequence 
v = ±2, ±6, ±10, • • • due to the Dirac fermion nature of 
the particles and the initial SU(4) symmetry [HQ. When 
disorder becomes weaker other more fragile QH states in- 
duced by interactions emerge. For instance spontaneous 
symmetry breaking gives rise to the states with SU(2) 
symmetry at v = 0, ±4, ±8, • • ■ as well as the fully broken 
symmetry states at v = ±1, ±3 ; • • ■ Also fractional 

quantum Hall states appear [9Ml3l|. 

A key requirement to observe these fragile states in 
transport has so far been the fabrication of better qual- 
ity samples. This has been accomplished by pla cing 
graphene on BN [HI, EH, or by suspending graphene 11, 
Il2| . In the quest for observing fragile or novel incom- 
pressible states, probing a smaller area may be an al- 
ternative to circumvent the challenges of producing even 
higher mobility samples, since the sample may be much 
cleaner on the nanometer scale IB"18|- However, to ob- 
tain such microscopic information we are usually forced 
to resort to sophisticated local probe methods which pre- 
vent the blurring by disorder induced averaging. Indeed, 
transport studies normally offer an averaged view of the 
sample only and this averaging masks the interaction in- 
duced physics. 

That said, we note that graphene magnetotransport 
traces are frequently cluttered with fluctuations. Here we 
demonstrate that these fluctuations reflect processes that 
occur very locally. They are a manifestation of charge lo- 
calization and are more pronounced in the transconduc- 
tance. A systematic study allows observing higher or- 
der fractional quantum Hall and broken symmetry states 
even if these do not show up as quantized states in the 
Hall or longitudinal resistance traces. They make phe- 



nomena on the nanometer scale visible in a macrocopic 
transport experiment despite significant disorder. 

Our studies were carried out on both supported and 
suspended exfoliated graphene in either a two-terminal 
or four-terminal configuration. The representative data 
discussed here stem from a suspended two-terminal de- 
vice shown in Fig. 1. Fabrication procedures followed 
closely those previously reported [11 The device 

had a width of 3 (im and was about 0.8 /im long. A dc 
source-drain voltage was applied. The density was tuned 
with the doped Si substrate separated from the graphene 
device by a 150 nm thick SiC>2 layer and an equally thick 
air gap. Charge neutrality was reached at a gate voltage 
Vb g of about 1.2 V. The density inhomogeneity, And ~ 
6 x 10 10 cm~ 2 , was determined from the full-width at 
half- maximum of the field effect resistance (not shown). 

Figure 1(a) plots the two-terminal conductance in 
units of e 2 /h as a function of the density for perpendic- 
ular fields from 1 T up to 15 T in 1 T steps. There are 
no surprises here. Plateaus or precursors of a plateau are 
observed for the integer states at filling v — ±2, ±6, • ■ • as 
expected for monolayer graphene. Also broken symmetry 
states at v — and v = ±1 can be discerned. They have 
been attributed to the lifting of the spin and/or valley 
degeneracy of the Landau levels. Other broken symme- 
try states are missing and also fractional states remain 
hidden. Figure 1(b) displays the single trace recorded 
at 15 T for closer inspection. It is not your textbook 
quantum Hall curve and looks noisy. We are not deal- 
ing with noise though. These fluctuations are amplified 
in the transconductance <? m . When applying an addi- 
tional ac modulation 5Vh g to the backgate, an oscillating 
source-drain current SIds is induced and g m — SId s /SVh g - 
The frequency of (5Vb g was 433 Hz and its root-mean- 
square amplitude was 10 mV. This corresponds to a den- 
sity change of 5.8 x 10 8 cm~ 2 or the average addition of a 
single electron to an area with 0.5 /Ltm diameter. In the 
transconductance, the background has vanished and only 
the fluctuations are left over in Fig. 1(c). This greatly 
simplifies a systematic study of these fluctuations. 
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Repeating this experiment for different B-values yields 
the color rendition of Fig. 2(a). The rich set of features 
can be classified into several groups of parallel lines. For 
instance at the center, we find vertical lines parallel to the 
v = 0-line. A bundle of lines parallel to v = —2 and v = 2 
appears left and right. At the bottom, lines with a slope 
identical to the v — ±6-line can be discovered. These 
lines parallel to filling factor slopes at which the system 
is expected to condense in an incompressible state are 
reminiscent of the spikes that have been observed in lo- 
cal compressibility measurements using a single electron 
transistor on graphene [l9[ as well as of earlier local 
compressibility data on GaAs two-dimensional electron 
systems [201 ] . They were attributed to non-linear screen- 
ing and charge localization in a landscape of compressible 
dots separated by an incompressible sea, which forms as 
we approach complete filling of a Landau level. Since this 
localization picture is crucial for understanding Fig. 2(a), 
it is described here. 

A widespread view of localization in the quantum Hall 
regime is based on the association of localized states with 
equipotential lines which either enclose potential hills or 
are trapped inside a valley of the disorder potential. Plot- 
ting the position in the density versus B-field plane at 
which such localized states get filled, would generally re- 
sult in curved, non-monotonic lines that may cross as 
they change their shape and move up or down in the 
disorder landscape when changing B to maintain a fixed 
number of enclosed flux quanta. None of these properties 
are compatible with the observed parallel lines. This pic- 
ture only holds in the strong disorder limit when single 
particle physics prevails [20j |. 

Localization proceeds in a different manner as illus- 
trated in Fig. 3(a). Screening and Coulomb blockade 
physics are key ingredients. Charge carriers will redis- 
tribute and generate a spatially dependent density pro- 
file in an attempt to flatten the bare disorder potential. 
As long as the local density of states is not exhausted, 
they are able to accomplish this task. As a Landau level 
approaches complete filling in some regions the required 
density to flatten the bare disorder potential locally ex- 
ceeds the level degeneracy. Hence, the bare disorder po- 
tential can not be screened away. Here a potential barrier 
emerges. Eventually, compressible regions or dots are 
enclosed by an incompressible lake, Coulomb blockade 
physics sets in. The discrete nature of charge becomes 
relevant. Elementary charges can be added to the re- 
maining compressible dots only one at a time when the 
overall density is raised sufficiently. Each charging event 
corresponds to filling a localized state. As the average 
density is increased further, more and more charges are 
added to the dots. The same landscape of compressible 
areas appears and the same charging physics recurs at a 
higher B, but at the same density deviation from com- 
plete filling. The filling of specific localized states there- 
fore evolves along lines parallel to each other and to the 



corresponding integer filling factor in the (n, B)-plane as 
schematically drawn in Fig. 3(a). The same holds for 
fractional states but with fractional charges . 

It is natural to attribute the parallel lines in the 
transconductance to charge localization physics in this 
network of dots that forms when charge carriers are lo- 
cally unable to screen the bare disorder. Current may 
flow from the source to drain through this network as 
well as extended edge states (inset Fig. 3(c)). During a 
backgate modulation cycle, the Coulomb blockade may 
be lifted for some dots or activated for others and trans- 
port channels are turned on, off or altered. This results in 
an oscillating source-drain current. At a different B, the 
same transport channels would appear and disappear at 
the same deviation from complete filling. Hence, a plot 
of the transconductance in the (n, £?)-plane should be 
crowded by sets of parallel lines running along filling fac- 
tors at which incompressible ground states form much the 
same way as in local compressibility measurements [l9j . 
This is indeed confirmed in the experiment of Fig. 2(a). 
Weak line features were previously also observed in GaAs 
2D electron systems 22| and graphene 23|, l24j when 
just measuring the resistance. They were connected to 
single particle integer quantum Hall states only. The 
transconductance here reveals the much more fragile bro- 
ken symmetry states and fractional quantum Hall states 
in graphene. For instance, lines parallel to v = 0, v = 1 
and v = — 1/3 are easily discovered without further data 
processing in Fig. 2(a). Additional incompressible states 
can be identified by using a numerical analysis. It is 
based on the calculation of the correlation function C(v) 
= J2(DkiD pq 5 v ) for the discrete data set. The calcula- 
tion is illustrated in Fig. 3(b). The data form a matrix 
and Dij corresponds to the data point at density and 
field Bj. The product of any two data points D^i and 
D pq contributes to the sum in C(v) provided they fall 
within a band in the (n, J3)-plane bordered by two lines 
with a slope corresponding to filling v and with a width 
of two pixels. In this case 5 V equals one. Otherwise, 5 V 
equals zero and this pair of points does not add to the 
sum. C(y) is evaluated as a function of v and the result 
is plotted in a polar diagram as shown in Fig. 2(b) — (c). 
The angle corresponds to the filling. To reduce the com- 
putational effort the calculation of C(y) is restricted to 
a limited rectangular window in the (n, £?)-plane. The 
center of the analysis window is placed in areas where in- 
teresting features are expected. Key locations have been 
marked in Fig. 2(a). C{y) for different analysis windows 
are overlaid, but colored differently. In Fig. 2(b) maxima 
of C(y) appear at v = —2, —4, —6 and —10. Note that the 
v = —4 quantum Hall state, attributed to a symmetry 
reduction from SU(4) to SU(2), did not show up in the 
conductance up to 15 T. In the color rendition lines for 
other broken-symmetry states at v — and v — ±1 are 
seen down to ~0.5 T and ~1 T, respectively while in the 
two-terminal conductance these broken-symmetry states 
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give rise to plateaus or minima only above 4 T. Fractional 
quantum Hall states at filling factors v = —1/3, —2/3 are 
clearly seen in the C(y) spectra of Fig. 2(c). They are 
the prominent fractional states observed in the cleanest 
of all graphene flakes [III E3- But a lso the higher or- 
der v = —2/5 fractional state is present in Fig. 2(b) as 
predicted in Ref. 0, [l(J and as recently reported in local 
compressibility studies [Hj]. A small feature is seen near 
v = 1/2, but more work is needed to firmly establish 
its significance. None of the observed fractional quan- 
tum Hall states appeared in the conductance. The ap- 
pearance of parallel lines in the transconductance serves 
as indirect evidence for screening governed localization 
and the emergence of a dot network. In graphene more 
direct evidence can be obtained. Under normal quan- 
tum Hall conditions the transport through this network 
is short-circuited by the extended edge channels. They 
prevent a measurement on only this network. However, 
graphene exhibits a true insulating quantum Hall state 
near v = due to valley symmetry breaking of the zero 
energy mode [5j-[8|]. In this regime, the chemical poten- 
tial does not cross any of the valley split Landau levels 
even at the edge where they bend up or down. An inset 
of Fig. 3(c) displays the level diagram. Source-to-drain 
current flow then can only proceed by charge tunneling 
through the dot network. The differential conductance, 
dlds/dVds, as a function of Vbg is plotted at the top of 
Fig. 3(c). At the border of the insulating regime, trans- 
port features are present. Deep into the v — state, the 
system is too incompressible and no current can flow. 
The bottom panel of Fig. 3(c) displays the differential 
conductance as a function of Vb g and Vd s . It resembles a 
stability diagram and exhibits a series of Coulomb di- 
amonds. This confirms the existence of a network of 
dots. From the diamonds it is possible to extract typ- 
ical dot sizes. They range from 160 to 400 nm and re- 
flect the disorder length scale. The conductivity near 
the charge neutral point at zero field is approximately 
3.5e 2 //i. From this value, we can deduce a mean free path 
^mfp ~ 350 nm 26], which is consistent with the length 
scale estimated from the Coulomb diamonds. The irreg- 
ularity of the Coulomb blockade oscillations in dl^jdV^ 
is attributed to the network character or to the irregu- 
lar shape of the QDs for which chaotic resonances have 



been predicted [27]]. We note that signatures of such a 



complex network of QDs in graphene have recently also 
been reported in scanning tunneling microscopy [l8j. 

The transconductance fluctuations in Fig. 2(a) allow 
the extraction of a lower bound for the density variation 
across the sample. When filling localized states produces 
features parallel to a nearby quantum Hall filling factor 
2/j, it confirms a posteriori that the landscape of dots is 
independent of the B-field and the assumption of a fixed 
bare disorder potential is correct. It also implies that 
the localized states associated with quantum Hall state 
Vi all get filled within a fixed density band whose width 



Arid is determined by the difference between the density 
extrema n max and n m ; n of the spatial density variation 
that emerges to screen the bare disorder potential. In 
the large B- field limit, these bands of localized states 
are well separated and their width immediately yields 
the maximum density variation across the sample. Here, 
charging lines running parallel to adjacent quantum Hall 
fillings and z^l still overlap. For instance, features 
parallel to fn = ~2 are seen beyond the z/l = — 1-slope. 
Hence, the density variation exceeds B/<pQ-(ua — z/l) even 
up to 15 T. This case is illustrated in Fig. 3(a). This sets 
the lower bound for And at 3.6 x 10 n cm~ 2 . This lower 
limit is much larger compared with the disorder estimate 
from the field effect curve (6 x 10 10 cm~ 2 ). The latter 
presumably reflects the density variance rather than the 
difference between the extrema. 

We conclude that transconductance fluctuations offer 
a window to fragile incompressible quantum Hall states, 
because they reflect charging effects on the nanometer 
scale due to the appearance of the gap in the energy 
spectrum even if the disorder is so large that the quanti- 
zation signatures of these Hall states are entirely missing 
in standard magnetotransport. The method described 
here is an attractive alternative in the quest for observ- 
ing unconventional quantum Hall states. 
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sample preparation. 



K. S. Novoselovet al., Nature 438, 197 (2005). 
Y. Zhang et al, Nature 438, 201 (2005). 
Y. Zhang et al. Phys. Rev. Lett. 96, 136806 (2006). 
M. O. Goerbig, R. Moessner, and B. Dougot, Phys. Rev. 

B. 74, 161407 (2006). 

K. Yang, Solid State Comm. 143, 27 (2007). 
V.P. Gusynin, V.A. Miransky, S.G. Sharapov, and LA. 
Shovkovy, Phys. Rev. B, 77, 205409 (2008). 
J.G. Checkelsky, L. Li, and N.P. Ong, Phys. Rev. Lett. 
100, 206801 (2008). 

J. Jung and A.H. MacDonald, Phys Rev. B, 80, 235417 
(2009). 

C. Toke and J.K. Jain, Phys. Rev. B, 75, 245440 (2007). 
Z. Papic, M.O. Goerbig, and N. Regnault, Solid State 
Comm. 149, 1056 (2009). 

X. Du, I. Skachko, F. Duerr, A. Luican, and E.Y. Andrei, 
Nature 462, 192 (2009). 

K.I. Bolotin, F. Ghahari, M. D. Shulman, H.L. Stormer, 
and P. Kim, Nature 462, 196 (2009). 
CR. Dean et al, Nature Phys. 7, 693 (2011). 
CR. Dean et al., Nature Nanotechnol. 5, 722 (2010). 
G. Li and E.Y. Andrei, Nature phys. 3, 623 (2007). 

D. L. Miller, K.D. Kubista, GM. Rutter, M. Ruan, W. 
A. de Heer, P.N. First, and J. A. Stroscio, Science 324, 
924 (2009). 

[17] Y.J. Song et al, Nature 467, 185 (2010). 



[1 
[2 
[3 
[4; 

[5] 
[6] 

[7 

[8 

[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 



[18] S. Jung et at, Nature Phys. 7, 245 (2011). 

[19] J. Martin, N. Akerman, G. Ulbricht, T. Lohmann, K. von 

Klitzing, J.H. Smet, and A. Yacoby, Nature Phys. 5, 669 

(2009). 

[20] S. Hani, J. Martin, E. Teitelbaum, J.H. Smet, D. Mahalu, 

V. Umansky, and A. Yacoby, Nature, 427, 328 (2004). 
[21] J. Martin et at, Science 305, 980 (2004). 
[22] D.H. Cobden, C.H.W. Barnes, and C.J.B. Ford, Phys. 

Rev. Lett. 82, 4695 (1999). 
[23] S. Branchaud, A. Kam, P. Zawadzki, F. M. Peeters, and 

A. S. Sachrajda, Phys. Rev. B, 81, 121406 (2010). 
[24] J. Velasco, Jr., G. Liu, L. Jing, P. Kratz, H. Zhang, W. 

Bao, M. Bockrath, and C.N. Lau, Phys. Rev. B, 81, 

121407 (2010). 

[25] B.E. Feldman, B. Krauss, J.H. Smet, and A. Yacoby, 

larXiv:1201. 51281 
[26] E.R. Mucciolo and C.H. Lewenkopf, J. Phys.: Condens. 

Matter, 22, 273201 (2010). 
[27] J.H. Bardarson, M. Titov, and P.W. Brouwer, Phys. Rev. 

Lett. 102, 226803 (2009). 



5 




FIG. 1: (a) Two-terminal conductance as a function of density at fields from 1 T to 15 T in f T steps, (b) — (c) Comparison 
of conductance (b) and transconductance (c) curves measured at B = 15 T. Both measurements were done with a dc-voltage 
bias of Vds = 500 fiV. Insets show the schematics of the measurements. 




FIG. 2: (a) Color rendition of the transconductance in the (n, B)— plane, (b) — (c) Correlation function C obtained by analyzing 
the data in the windows shown in (a). The details of this analysis method are shown in Fig. 3. 
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FIG. 3: (a) Schematic illustration of the formation of QDs 
and the development of the compressible spikes in the (n — B)- 
plane. (b) Schematic illustration of the numerical analysis to 
obtain the spectra shown in Fig. 2. (c) Differential conduc- 
tance and stability diagram measured at 1.4 K and 14 T. 
Inset shows an energy diagram along the sample on the left 
and a schematic of the network of of compressible QDs with 
extended edge states on the right. 



